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Abstract 

Newtonian Cosmology is commonly used in astrophysical problems, be- 
cause of its obvious simplicity when compared with general relativity. However 
it has inherent difficulties, the most obvious of which is the non-existence of 
a well-posed initial value problem. In this paper we investigate how far these 
problems are met by using the post-Newtonian approximation in cosmology. 
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1 Introduction 



In cosmology the use of Newtonian hydrodynamics is frequently favoured over the 
more correct theory of general relativity ^, ^. While the main reason for this is 
its obviously greater simplicity, there are aspects of the Newtonian approximation 
which have distinct disadvantages. In the first instance, it is not always understood 
what actually constitutes the Newtonian approximation to general relativity. Al- 
though generally assumed to be a specialization of the linear approximation, this is 
not strictly true for the Newtonian hydrodynamic equations are not the linearized 
approximation of the Bianchi identities, since the term involving the gradient of 
the potential can only be obtained by considering a higher (non-linear) order of 
approximation. 

Secondly, and not totally unrelated to this difficulty, is the fact that the equations 
of Newtonian hydrodynamics with gravitation do not have a well-posed initial value 
problem. Without a well-posed Cauchy problem the future of initial fluctuations 
of density and velocity flelds cannot be uniquely determined from the Newtonian 
equations, yet a major reason for using Newtonian theory in cosmology is that its 
perturbation theory, both linearized and exact, is considerably easier than that of 
general relativity 0, . What value are we to place then on the results of Newtonian 
perturbation theory and its consequences for galaxy formation? 

On the other hand it is well known that general relativity has a well-posed 
Cauchy problem for perfect fluids with an equation of state of the form P = f{p) 
[0, 1^ . What then goes wrong in taking the appropriate limit to Newtonian theory? 
The principal aim of this paper is to write down a higher order approximation of 
general relativity which results in a closed system, including Bianchi identites, and 
has a well-posed Cauchy problem. Although higher order approximations in general 
form an endless sequence of equations |]T0|, |1T[, called successively post-Newtonian, 



post-post Newtonian etc., we will show that possible to make the cut-off at the 
post-Newtonian level in such a way that these goals are achieved. 

In Section 2 we outline the standard theory of Newtonian Cosmology and de- 
scribe in detail the difficulties stated above. We then derive in Section 3 a standard 
sequence of approximations to general relativity. It will be seen that wherever one 
curtails this sequence the Cauchy problem remains ill-posed. In Section 4 the equa- 
tions of the full post-Newtonian approximation are then cast in a form which gives 
a closed system having a well-posed Cauchy problem. In Section 5 we discuss future 
plans for this theory. 

2 Difficulties with Newtonian Cosmology 

Following the traditional approach of Heckmann and Schiicking ^, a Newtonian 
cosmology is defined to be a three dimensional Euclidean space and a universal time 
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parameter t, with a perfect fluid matter source for a scalar gravitational potential 
field. Defined on the spacetime are three scalar functions, the density p(r, t), pressure 
P{r,t) and gravitational potential (j){f,t), and a velocity field vector v{f,t). These 
are related by the standard equations of fluid dynamics, namely the equation of 
continuity, Euler's equation, and Poisson's equation: 



where throughout this paper a dot will refer to the partial derivative d/dt. 

The system of equations (|l])-(|D together with an equation of state P = f{p) is 
not well-posed however, as there is no proper initial value problem. Suppose one is 
given initial values po{r) = p{f, to) and Vo{r) = v{f, to), then clearly Po{r) is obtained 
from the equation of state, and the gravitational potential 0o(^ can be found by 
solving equation (^. Note however unlike standard Newtonian mechanics where the 
gravitational potential and its derivatives are always assumed to vanish asymptoti- 
cally at large distances, in Newtonian cosmology no obvious boundary conditions at 
spatial infinity suggest themselves since the density is no longer assumed to vanish 
at infinity For this reason, the potential is only known up to an arbitrary solution 
of Laplace's equation V^V = (harmonic function). While the time derivatives po 
and vq are determined from equations (|T]) and (|^), the time derivative of must 
also satisfy a Poisson equation. 



Again there arises an arbitrary harmonic function ipi{f). This process continues 
for every higher time derivative of 0, giving rise to an infinite number boundary 
condition problems. The system of equations is not well-posed since t = const is 
a characteristic surface of the system and it is not permitted to set initial data on 
such a surface 0. 

At first sight this is a peculiar phenomenon, since Newton's theory is a limiting 
approximation of general relativity, yet Einstein's equations are well known to have 
a well-posed Cauchy problem. Some structural information is clearly being lost in 
taking this approximation. 

The usual Newtonian theory is obtained from general relativity by linearizing 
gravity and taking small characteristic velocities. In the weak field limit the metric 
has small deviations e^jy from flat space and takes the form 



ij + v-Vv = -V0 VP, 

p 



p + V- Vp + pV ■{7=0, 



(1) 
(2) 



= AnGp. 



(3) 



= AirGp. 
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As usual Tjfj^i, = diag(— 1, 1, 1, 1), Greek indices run from to 3 and Latin indices 
from 1 to 3. 

Since g^^gx^ = 5^, it follows that 

g^"" = 7]^"' -e^''' + 0{e^), where £'^'^=77^^%^. 

Under an infinitesimal coordinate transformation 

with corresponding gauge transformation 

it is possible to choose coordinates for which the harmonic gauge condition holds, 
or equivalently, 

= 0, (4) 

where ^ 

In the harmonic gauge, the linearized Einstein equations read 

G^y = -^y^tiu'"^^ = kT^i/ where k = 8ttGc~'^ (6) 
For a perfect fluid the energy-stress tensor takes the form 

T^, = {pc^ + P)U,U, + Pg,,. (7) 
where, in the Newtonian approximation, the 4-velocity has components 

f/. ^ (-1, -) +0(/32), where /3 = < 1. 

c 

The pressure P will be assumed to be of order 0{l3'^)pc^ since it is approximately 
the kinetic energy density, whence to highest order in each component 

Tm = pc^, = Tio = -pViC, Tij = pViVj + PSij. 

Consistency in Einstein's equations (P) can only be achieved in the limit /3 — > 
if it is assumed that e and (3 are the same order of magnitude, and 



where is the Newtonian gravitational potential, assumed to be of magnitude 
0(/3^c^). This is equivalent to assuming 

ip^, = diag(-40c-2, 0, 0, 0) + 0{(3^). 

Finally, if for all physical quantities, derivatives d/dx^ are assumed to be of or- 
der f3d/dx^, then the (0,0) component of Einstein's equations (|^) gives Poisson's 
equation (^. 

At this point we meet the first difficulty. The remaining Newtonian equations 
(|l]) and (^ should be obtained from the Bianchi identities. However the linearized 
Bianchi identities are not consistent with these equations, since on using the har- 
monic gauge condition (^) it follows immediately that T^'^^,^ = = 0. Sub- 
stituting the above expression for the energy-stress tensor then results in an Euler 
equation without the term p0^j. The usual resolution of this problem is to include the 
term FqqT^" ^ —|pc^ eoo,i arising in the full (non-linear) Bianchi identites T^'^.^y = 0. 

This mixture of linearized and non-linear approximations to general relativity is 
clearly unsatisfactory. It will be shown in Section 3 that by including the highest 
orders of approximation for the other 9 Einstein equations, which have so far been 
ignored, a version of Newtonian theory can be arrived at which yields the correct 
Bianchi identities. While less familiar, this version of Newtonian theory contains 
precisely the same information that is present in equations (|ll)-(3)- It therefore 
resolves non of the difficulties discussed above. 

We conclude this section with a brief discussion of the Cauchy problem in general 
relativity, and why it is that the Newtonian limit loses the well-posed character of 
Einstein's theory. Following the treatment given by Synge ^ with minor modifica- 
tions, a suitable initial data set given at x° = is g^i,, gij^o, p and Ui. The pressure is 
assumed given by an equation of state P = f{p). In order to obtain a unique solution 
of G^u = i^T^u it is necessary to subject this initial data to constraint equations 

= (8) 

and the evolution equations 

g'^K. = 0, (9) 

= 0, (10) 
Rij = K{Ti,-^Tgij) where T = T^. (11) 

Equation (^ is the full harmonic coordinate gauge condition and provides evolution 
equations for (/q^, while the Bianchi identities (|^) are evolution equations for p and 
Ui which ensure that the constraint equations are maintained at later times. The 
evolution of the metric components gij are determined by Eqn ([TI|) which has the 
structure 

~ ~2g gijfio + ID, 
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where ID refers to terms expressible entirely in terms of initial data quantities. 

When the Newtonian approximation is taken, ^f^^, ^ rj^y — 2(f)c~'^6f^u, the only 
remaining equations at highest order 0(/5^) are as follows: the constraint equations 
reduce to the single Poisson's equation (H), the full Bianchi identities approximate 
to equations (|ID and (H) while the Rij evolution equation reduces again to Poisson's 
equation. Note how this latter evolution equation has become identical with the 
constraint equation in the Newtonian limit. It is also worth observing that the full 
Bianchi identities are necessary to give the correct Euler equation. 



3 The Newtonian Approximation of a General 
Relativitistic Perfect Fluid 

Following a schema similar to that adopted by Weinberg |10], we adopt units in which 
the typical velocity has magnitude 1, i.e. (3 ~ c~^, and assume a one parameter family 
of metrics g^^^x^^ c) for which there is a system of coordinates (a;°, x*) in which the 
components have the following asymptotic behaviour as c — > oo: 

= -1 - 20c~2 - 2ac-^ - 2a' c-^ - 2a" c'^ , 

goi = CiC-' + C^^' + Cc-' , (12) 

9ij = ^ij - '^(pdijC'^ + aijC~^ + a'ijC'^ + oi'lf'^ 



It will also be useful to expand the quantity ip^u defined in 

ipoo = -4(f)c-^ + ec-^ + e'c-^ + e"c-^ 



whence 



1 1 

oiij = <Pij + ]j^kj{.d - (pkk), oi'ij = (p'ij + - <Pkk), ■■■■ 
The Harmonic Gauge condition @ gives a series of equations in successive powers 

of C-2, 

^ = (13) 

Ci = <Pij,j, (14) 

= C'„ (15) 

a = 0:,,,, (16) 
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Replacing the quantities which define the perfect fiuid by the expansions 



p — > p + p'c'^ + p"c-^ + ... 



P P + P'c'^ + P"c'^ + ..., 

\Uo\i-l, ViC-^+v[c-^ + v'^c-'' + ...), 

it follows from U^U^ = —1 that 

Uo^ = l + (20 + v^)c-^ + {2a + 6(f)v^ + + 2CiVi + 2viv[)c-^ + ■■■. 

The energy-momentum tensor is defined (|^) as before, with an equation of state 
P = f{p), which implies for higher pressure terms, 

P' = f\p)p, P" = f\p)p" + \f"{p){p")\ ••• (17) 

where primes on the function / refers to its derivatives. Expanding the Ricci tensor 
as a power series in c~^, 

and substituting in the Einstein field equations 
we find in the harmonic gauge, to order c~^. 



where 
and 



= 47rGp, (18) 

Ci,kk = IGTiGpVi, (19) 

4>ij,kk = -lQ7rG{pViVj - 6ijP) + Aij, (20) 

^%oo: = \i-e,kk + A), (21) 



A,j = 800,,j + 40,,0 J - 6ij (60,fc0,fc + 32nGp(P) , (22) 



A = 60,,0,i- 167rG(pt;2 + 4p0-p')- (23) 
The first three equations of this set (p!8D -(^0D together with the harmonic gauge 



conditions (0) and (p^ ) constitute a reformulation of Newtonian cosmology since 
the time derivate of ([T^) and d/dx^ of ([T9| ) together give rise to the equation of 
continuity (|1]), while d/dt of (|19|) and d/dx^ of ( pOD give rise to Euler's equation (0) 
with the correct gravitational term. 

This version of Newtonian cosmology seems a little strange, since it introduces a 
new vector field Q and tensor field It should be realized however that these two 
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fields play an entirely subsidiary role in that the equations of Newtonian cosmology 
(HD 5 (i) 3) precisely the integrability conditions for (|l^) , ([l^) , (|19|) and (pOD . 
The fields Q and (pij play no further part in the theory. In similar vein, although 
our version of Newtonian theory is a 4th order approximation to general relativity, 
the remaining 4th order equation (^TJ) is also entirely subsidiary in nature, since it 
merely serves to define the next order approximation to density, p', in terms of and 
arbitrarily specified quantity 9. 

Although equations (|13D and (|14D give rise to a time evolution equation for the 
gravitational potential, 

this in no way helps with the well-posedeness problem since in order to obtain higher 
time derivatives of (p it will be necessary to solve an infinite sequence of complicated 
Poisson-like equations for higher derivatives of (pij which arise on taking successive 
time derivatives of eq. (pO|) . 



4 The Post-Newtonian Approximation 

Continuing the approximation of Einstein's equations to order c"^, results in the 
equations 



Ci 



Oij 



1 1 



Vij,kk + Bij, 

0[kk + c, 



(24) 
(25) 
(26) 



where 

B, 

B.;a 



pv[ + Vi{p' + P) + pViV - ^pCi 

2 idCk,kj ~l" CjCk,ki) Ck {Ci,jk ~l" Cj,ik) ~l" '^CkCk,ij Ck,iCk,j ~l" Ci,kCj,k 
-20_A. {(f)ki,j + <Pkj,i - 20ij- fc) - 1600,i0j + 20_i6'j + (f),i4>mm,j + 20j6',i 

+(p,j4>mm,i — 20 {(pkijk + 'PkjAk — 4>ij,kk — 4>mm,ij — O^ij) — 2(j)ki(t>Jk 

\2 



(27) 



— '^4>kj(p,ik + 20,jj(6' + (prnm) — ^ij jCk,iCk,i + \Cj,kCk,j + |(Cfc,A:)^ 
-CkCm,mk - ^4>,k4>ki,i + 40,fc0ii,fc - 1200,fc0,fe + ^,iO,i - 2(j)ki0,ik^ 

-0 (20fci,ifc - 2(l)u,kk - + SttG \2{p' + P)viVj + 2p{viv'- + Vjv'^) 



+2p(20 + v^)viVj + pcPij + 5ij {2P' + 2pcPv^ - ^pO - \cPp' - \ 
+ j(l>,kP,k + hp(l)kk 



(28) 
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and 



C = ^Ck,iCk,i + |(Cfc,fc)^ ~ ^Cj,kCk,j + '^4>,k4>ii,k — '^^4>(l> ,k4> ,k + 

— 2(j)ki6,ik — 4> {^4>ki,ik + '24>ii,kk + l^ii) + ^4>ki4>,ik 

.2 9, 



+87rG \2p" + |P0 - |0p' + 2p't;2 + 2Pt;^ - |0,fcp,fc 
+p (4t;^t;, + Qv^<P + 2v^ - \<Pii + ^6 - 80^ + 4^;^^) ] • (29) 

Equations (p5| ) and (p6| ) require a certain amount of juggling before they can be cast 
in form given here. While eq. (^) does provide an equation for 0^^, this inviting 
feature is counteracted by the fact that cfy'^^ is not arbitrary but must satisfy the 4th 
order constraint equation 



V'Bij + A. 



16nG 



{pviVjy+ 6ijP 



where all time derivatives on the right hand side can be reexpressed in terms of ini- 
tial data quantities (undotted quantities, but possibly involving spatial derivatives) 
via the various evolution equations such as the harmonic conditions and Bianchi 
identities. Without this constraint equations (|2T|) and (^) will not be consistent 
with each other. Consequently all that has been achieved is to push the system to a 
higher order of accuracy, but with no further resolution of the initial value problem. 
Reading off successive approximants by "peeling off" the higher powers of in 
Einstein's equations continues the same problem to higher and higher levels. At no 
stage is it possible to close off the system in a self-consistent way. 

There is however another method [|. Instead of reading off each successive power 
as a separate equation, suppose we "chop off" Einstein's equations @ at successive 
levels c~^, etc., retaining the entire equation instead of peeling off its component 
parts. The order theory is then the Newtonian approximation given in the 



previous section, while the order c 
wavelike equations 



theory replaces equations 



by the 



— c 



',kk 







2a 
C si,kk 



-A7rGpc''--{9,kk-A), 
-IQnGpViC^ + C kk + Bi 



(30) 

(31) 
(32) 



(pij - c (pij^kk = 4>'ijMk + Bij + c [167iG{pviVj + 6ijP) - Aij] , 

where the expressions for A, Aij, Bi and B^j have all been defined above. 

Using the harmonic condition (|T^ and (0) we obtain the following Bianchi 



identities from (|30D , (|3TD and (|3 



P + {.PVi),i + 



167rGc2 



kk 



C 



jjkk 



A-B 



J.J 



(33) 



^We are indebted to D. Hartley for this proposal. 
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p{Vi + VjVij + + Pi + (p + {pVj),j)Vi 

= + - - + c'A.^j) (34) 

These equations form a closed and well-posed system in the following sense: 

1. Set 10 arbitrary functions of space and time, (p'^if^t), p'{r,t) and f-(f, t). 

2. Rewrite eqs. (pOD and (^) as a set of 4 constraint equations by substituting 
= —\(t>ij,ij and Ci = which follow from the gauge conditions (|l^) and(|l^, 



kk 



WnGpVi + 



\kk 



„2 ( 



Vjk.jk 

^kk + 4 



A), 
-B, 



(35) 
(36) 



3. Finally there are a total of 18 evolution equations, consisting of ([T3|), (p^, 
(H), (P, (H, (H), and (II- 

The system is well-posed, for given initial data (pij{f, 0), (pij{r, 0), 9{f, 0), CK?^, 0), 
p(r, 0), fj(r, 0) all arbitrary functions of r, and 0(r, 0), 0(^)0) subject to the con- 
straint equations ( pS]) and (|36|), then the evolution equations determine a unique 
space-time dependence for cf), Q, (pij, 9, p and fj. What makes the system closed 
and self-consistent is the fact that ( p3D and the , i derivative of (|36| ) imply d/dt 
of the first constraint equation (^), while (Q) and ,i of eq. (|32D imply the time 
derivative of (P^). Hence the constraint equations are carried forward in time as a 



consequence of the evolution equations and will automatically be true at later times 
if they hold at t = 0. It should be remarked finally that eq. (^) is redundant in that 
it can be regarded as merely acting to define p" in terms of 6' and other quantities. 

Although every choice of functions 0^^, p' and v'^ leads to a viable post-Newtonian 
cosmological theory, it is obviously simplest to set them all to zero. It also follows 
then from ( |T7|) that P' = 0, while consistency with the harmonic conditions ([15|) 
and ( [T6| ) is most simply maintained by setting ^ = and Q = 0. With these 
simplifications the Bianchi identities ( ^31) and (|3^) read 



P 1 + 



40' 



{pv 



1 ^ \ 1 
1 + -7 +-T 



p (2vjVj + 2vjVkVk,j + ^Cj,j 



(20,iO,ij ~ 2Ci0jii - 3Ci,i0jj) 



(37) 



and 



p(,Vi + VijVj + 0,i) + 



1 



-{A + Bjj)v,+ 



Ba — B, 



20(A j + 0jyt,fci) - 0i(y4 + (t)-jk.jk) 



(38) 
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5 Conclusion 



In this paper we have shown that standard Newtonian cosmological theory is inad- 
equate; the Bianchi identities are not obtainable from the field equations and there 
is no well-posed initial value problem. By keeping higher order terms it was shown 
that Newtonian theory can be reformulated as a new theory where the Bianchi iden- 
tities are consistent with the field equations. However this new Newtonian theory, 
although consistent, still does not have a well-posed initial value problem, and it is 
necessary to go to the post-Newtonian level in order to achieve a physically viable 
cosmological theory. 

The usual Newtonian gravity leads to models which have significant dissimilar- 
ities to the corresponding general relativistic models ||^. In particular the theory 
of anisotropic homogeneous cosmologies is quite different in the two theories. In 
Newtonian cosmology these models can, in the case of dust P = 0, all be classified 
by giving 5 arbitrary functions of time (the components of shear), while in general 
relativity it is necessary to classify models into the well-known nine Bianchi types. 
It will be interesting to see if this new post-Newtonian theory can produce a theory 
of homogeneous models more in line with the general relativistic scheme. 
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